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Abstract 

We provide an analytical calculation at the two-loop level of the real non-abelian 
contribution to the leading {1/Q) correction to the mean value of the C param- 
eter in e + e~ annihilation, to complement the existing calculation of the abelian 
contribution; and we compare the result with the numerical 'Milan factor' ob- 
tained using the soft approximation. We find agreement with the previous 
results. The use of the pinch technique to separate the various contributions 
yields insights into the structure of renormalon-type diagrams in a non-abelian 
theory. 
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1 Introduction 



The study of event shape variables in e + e~ annihilation and in deep inelastic scattering 
continues to provide important insights into the structure of QCD. On the theoretical side, 
the state-of-the-art perturbative predictions comprise exact next-to-leading order calcula- 
tions and (for the distributions) a resummation of large logarithms. In addition to the 
predictions of perturbation theory, large non-perturbative contributions must also be in- 
cluded. These often take the form of power corrections, which are terms behaving as 1/ Q n 
where Q is the hard scale of the process and n is some power. 

It has long been known that totally inclusive QCD observables, including for example 
the total e + e~ cross-section, possess operator product expansions that may be applied to 
the calculation of non-perturbative contributions However, no such expansion is avail- 
able for event shape variables, so some physically reasonable (and experimentally testable) 
assumptions are required before progress can be made. Either renormalon methods are 
employed (for a review see 0) with the assumption of ultra-violet dominance ||, or else 
assumptions are made concerning the behaviour of the non-perturbative strong coupling in 
the infra-red, leading to the dispersive approach Q . In practice the two methods give con- 
sistent results, but in the calculations that follow the dispersive approach will be adopted. 

We consider the resummation of loop insertions into a gluon propagator, which gives rise 
to a QCD effective charge. (Strictly speaking, this is valid only for fermionic loop insertions, 
but gluon loops may also be included using 'naive non-abelianisation' @, §J, or using the 
1-loop pinch technique J?) -[II).) This perturbatively-defmed effective charge diverges in the 



infra-red near the Landau pole, leading to an ill-defined perturbative prediction. But if it 
is postulated that the full non-perturbative coupling does not have this pole, but rather is 
well-defined and analytic throughout the complex plane except for a branch cut along the 
negative real axis, then well-defined predictions can be made. A dispersive representation 
for this coupling is adopted, in which the dispersive variable plays a role formally similar 
to that of a gluon mass. The non-perturbative contribution is defined to be the difference 
between the full prediction and that already included in the fixed-order perturbative result: 
for event shape means the leading correction is found to be proportional to A\/Q, where 

C F r , 2 s , 2, 



Ai = — / dfi nSa s (id ) (1.1) 

Z7T JO 

is the first half-integer moment of the non-perturbative component of the coupling, 5a s (fj, 2 ). 

It was however quickly pointed out (12| that the direct application of this method fails to 
take sufficiently into account the details of the gluon splitting, and that a more sophisticated 
application is required which treats this splitting correctly. Exact calculations in the abelian 
limit have been performed for the 1/Q corrections to the longitudinal cross section [|Tj|] and 
the mean value of the C-parameter |L4j] in e + e~ annihilation, while a more general analysis 
of the most commonly used variables, including non-abelian terms, was performed using 
the soft approximation in |L5|- ||17| . It is found that the effect of non-inclusiveness is simply 
to enhance the amplitude of the 1/Q corrections by a universal 'Milan factor', thus leaving 
the universality pattern unchanged. The extraction of the non-perturbative parameter A\ 
(or, more usually, the related quantity called a ) from a variety of event shape means and 
distributions is thus a stringent test of the model (as performed for example by [18|]). It is 
pleasing to find that universality of the coupling approximately holds. 

In the present paper we extend the exact e + e~ C-parameter calculation to include all 
non-abelian real emission terms, using the 1-loop pinch technique as a convenient tool 



1 



for the separation of contributions. Virtual contributions are not calculated here, but are 
included in the same manner as in [HJ. We recover the Milan factor of JK], which we 
present in closed form. 

The layout of the present paper is as follows: the relevant existing results are recalled 
in section [2|, while the calculation of the real non-abelian diagrams is presented in section 
[3]. A discussion follows in section |], in which we include the contributions from virtual 
diagrams and comment on the results. 



2 Recollection of existing results 



2.1 The C-parameter 

The C parameter |19| is defined by 



q n I 11 1 

^E^Fsin 2 % (2.1) 

where the sums run over all final-state particles, with p; being the 3-momentum of the ith 
particle in the centre-of-momentum frame of the collision (i.e. J2iPi = 0), and 0y is the 
angle between and pj. The hard scale is Q = J2i I Pi I- 
An equivalent frame-independent definition is 

C = 3 3 & ' Pj ^ (2.2) 

where the pi are the 4-momenta of the n final-state particles, and q is that of the decaying 
virtual photon, q = J2iPi- It is seen that C = for a pencil-like event, and C = 1 for 
a spherical event. In particular, C vanishes at Born level, so the first non-trivial order in 
perturbation theory is 0(a s ). 

At first order in a s , let p±, P2 and k be the respective momenta of the outgoing quark, 
antiquark and gluon. Then C takes the form 

C = 6 (1 ~ Xl)(1 ~ X2)(1 ~ X3) , (2.3) 

£1X2X3 

where = %Pi,2 • q/q 2 and X3 = 2k ■ q/q 2 are the energy fractions carried by the quark, 
antiquark and gluon respectively, and J2i x i = 2 by energy conservation. The leading-order 
result for the mean value of C is then |19| 



(C} = - [ daC= — a s (4vr 2 - 33) + 0(a 2 s ) . (2.4) 
a J 2n 

In order to calculate power corrections in the dispersive approach, we need to consider 
the gluon to be slightly off-shell. First we use the massive gluon scheme, which treats the 
gluon decay inclusively, and then we consider in more detail the splitting of the gluon into 
a quark-antiquark pair or into two gluons. 
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2.2 Massive gluon scheme 



Suppose we estimate the 1/Q power corrections to the mean value of C using the massive 
gluon prescription. Although this is known not to be a complete calculation, since it treats 
the gluon decay inclusively, it is a good first approximation and a necessary starting point 
for the full calculations below. 

If the outgoing gluon is slightly off-shell, with virtuality k 2 = eQ 2 , the expression for C 
becomes 



C = 3-3 



[I-X3 + e) 



x 2 



s) 2 (l- Xl - e f 



x x x 2 X3X1 
Then the expectation value of C is given to 0(a s 

1 



x 2 x 3 



4 



(2.5) 



by 



a 



daC 



167TiV r JO 



dxi 



1— xi- 



dx 2 W^W^C 



(2.6) 



where (—ie)(—ig)W afJ- is the tree level matrix element for the decay of the virtual photon 
with polarisation index /1 into a quark, an antiquark, and a gluon with polarisation index a, 
and we implicitly sum over fermion spins. Current conservation then implies that q^W^ = 
and k a W^ a = 0. 

In terms of the variables Xi we have 

"(xi + e) 2 + {x 2 + e) 2 e e 



W^ a W* a = 8N C C F 



and so from (12.61) we obtain 



1 - Xi)(l - x 2 ) 



[l-Xi) 



{C)r 



47r 2 - 33 - 127rv^+0(e) 



' mg 2tt 

Note that in the limit e — > we recover the first-order perturbative result 
sfe term gives rise to a leading power correction of 



5(C) 



mg 



12- 



Q 



(2.7) 

(2.8) 
), while the 

(2.9) 



where the non-perturbative parameter A\ is given by fll.l]) . 

In fact we obtain the same result whether we use the massive definition of C, ( ^.5| ), or 
the massless definition ( |2.3|) with the mass entering only via the phase space ( p.6|) . This is 
not always the case: in some cases the two conventions give different results, for example 



the longitudinal cross section <j£ in e + e annihilation [20|-|[22|| or the C parameter in DIS 
|23| . The Milan factor formulation is based on the massless shape variable with massive 



phase space, so this is the approach we must take when we wish to make comparisons. 



2.3 Abelian gluon splitting 

Suppose now the virtual gluon splits into a "secondary" quark-antiquark pair, with 
momenta ki j2 , as in figure [I]. The mean value of C now becomes 

2 T„r> r _ _ rl-e 



dk TjiUf 



J dU V s[k -> k h k 2 ] Tr^Ja^C . (2.10) 
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Here we have factorised the four-particle phase space into the massive-gluon phase space 
over X\ and X2, the Lorentz-invariant phase space for the gluon splitting, and an explicit 
integral over the gluon virtuality k 2 , which plays the role of dispersive variable. The 
running coupling a s (—k 2 ) is generated by making renormalon insertions into both gluon 
propagators. 




Figure 1: Diagram showing the virtual gluon splitting into a quark -antiquark pair. 



For C we take the full expression (|2.2| ) for four outgoing particles, which we decompose 
into four pieces: 

C = C (3) + C {p) + C {m) + C {s) . (2.11) 

The first contribution is just C^> = 3, the contribution from the first term of ( |2.2|) , which 
gives just three times the total cross-section. C^> contains the remaining "inclusive" terms, 
which are insensitive to the details of the gluon splitting and involve the momenta of the 
primary fermions only: 

(pi ■ q)(p2 • q) ' 

is the term involving the momenta of the secondary fermions only: 

r ( s) _ o • fc 2) 2 . . x 

° ~ V?)(Vg)' (2 ' 13) 

while is the sum of all mixed terms. Using the symmetries of the integral we may 
replace k 2 by k\ and p 2 by p\ such that 

C (m) = _ 12 (Pl^l) 2 , . 

We also make manifest the connection with the dispersive approach by writing for each 
of the contributions 

(C^U = %\% ^M-k^C^/Q 2 ) . (2.15) 

The factor (T R rif/3ii) \a s (—k 2 )\ 2 is none other than p s (k 2 ), the spectral function in the large 
rif (Abelian) limit. Thus we can compare with the 'naive' dispersive treatment, with the 
functions CW corresponding to the usual characteristic functions. 
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These were calculated in [14], where the total contribution to the C parameter from 



abelian gluon splitting was found to be: 



4 c i7r 3 

C ab (e) = 4vr 2 - 33 - + 0(e) . (2.16) 



Again we retrieve the standard result ( |2.4| ) in the limit e — > 0, but the coefficient of 
•y/e, giving the magnitude of the l/Q power correction, is modified by the full treatment of 
the gluon splitting. Indeed the exact calculation generates an enhancement of the naively- 
calculated result by a factor 157r 2 /128 = 1.157. 

3 C-parameter from e + e~ —> qqgg 

Now let us consider the splitting of the virtual gluon into two real gluons. In order 
to achieve gauge invariance we must include all tree-level diagrams that contribute to the 
process 7* — > qqgg: we cannot just take the loop insertion diagram shown in figure 0(a). 

For the purposes of calculation, some prescription must be chosen concerning the gauge- 
dependent contributions to the diagrams; this usually takes the form of a definite gauge 
choice. It is helpful here to use the one-loop pinch technique on the squared Feynman 
diagrams, in order that the cut loop insertion diagram 0(a) generates the non-abelian con- 
tribution to /3q, behaving analogously to the abelian contribution above, with the remaining 
diagrams giving additional contributions. Indeed if one recalls the result of inclusive inte- 



gration over soft gluon emissions given in [15 



1 ^r^^^^A-^^^^)' 



2! Jo Jo 2vr V w' W m 2 (k 2 ± + m 2 )\ ^ M m 4 

with m the invariant mass and k± the transverse momentum of the gluon pair, we see that, 
using the pinch technique, the contribution from figure |2|(a) generates the term proportional 
to j3o, while the diagrams like those of figure |](b, c) correspond to the additional logarithmic 
divergence. 

Note that, for the purposes of this calculation, the pinch technique is used simply as a 
means of reassigning gauge-dependent terms between diagrams, as explained in appendix 
[A], and no particular significance is attributed to any particular subset of terms. 

3.1 Non-abelian splitting (I): contribution to (3q 

The contribution to (C) from the cut loop insertion diagram, figure [|(a) using the pinch 
technique, is given (see appendix |A]) by: 

1 r dk 2 C a rl-« 1-3=1 -e 

(C) nal = - _^| as (_P)|2 / dxi / ^ dx2 w^W*/ / dLips[k - k u k 2 ] 

21 J K 07lly c JO Jl-xi-e J 

x {8(k 2 g af3 - k a k p ) + 2(k la - k 2a )(k w - k w )}C . (3.2) 



This is very similar to equation ( |2.10| ): apart from the symmetry factor 1/2! and the 



different colour factor, the only difference is in the coefficients of the two terms inside the 
curly brackets — recall that in ( |2.10| ) we had: 

Tr [7^17^2] = -^{k 2 g a f3 - k a kp) - 2(k la - k 2a )(k 1 ^ - k 2 p) ■ (3.3) 
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Pi 

Figure 2: Squared Feynman diagrams for the process 7* — > qqgg: (a) is the cut loop 
insertion diagram, with two internal gluon propagators of momentum k, (b) and (c) are 
examples of other diagrams, with one and zero such propagators respectively. 
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Thus the integration proceeds in precisely the same manner as in section |2.3| , yielding 
the result: 

<C W W = -%\ , (3-4) 
where the total characteristic function is now 

26l7T 3 

Cnai(e) = 4tt 2 - 33 - —^Ve + 0(e) . (3.5) 

lib 

For completeness the individual components C^i (e) are listed in appendix |B|. 

The factor (— 11(7^/12^) |a s (— /c 2 )| 2 is the 'spectral function' for rif = 0, which is nega- 
tive in this case. The running coupling a s (— k 2 ) is generated by making renormalon inser- 
tions into both gluon propagators, as it is possible to do using the one-loop pinch technique 
P"T|| . As before, the characteristic function reproduces the first-order perturbative result 
( [2.4D in the limit e — > 0, and gives a modification to the naive value of the 1/Q correction 
because of the non-inclusive nature of the observable. 

This contribution gives an enhancement of 877r 2 /704 = 1.220 over the 'massive gluon' 
calculation: the reason this is different from that found in the splitting g — > qq is because 
the splittings have different geometry, which in turn is due to the different coefficients in 
equations ( |3.2| ) and (|3 . 3|) . 



3.2 Non-abelian splitting (II): logarithmic divergence 



The diagrams that give rise to the logarithmic term in equation (3T) are those with 
only one or no internal gluon propagator, of which examples are given in figure @(b, c). 
These diagrams are considerably harder to evaluate than the cut loop insertion diagrams, 
for two reasons. Firstly, they are collinear divergent: an internal quark propagator can be 
made to go on-shell when one of the real gluons becomes collinear with either the quark or 
the antiquark. Thus we need a regulator, so we introduce a quark mass m 2 where 

m 2 JQ 2 = A < e (3.6) 

and we neglect terms that vanish as A — > 0. The divergent terms must however cancel in 
the sum of all real and virtual diagrams. 

Secondly, the squared matrix element does not as previously factorise into a tensor 
W m W^ independent of h\ and k 2 arid a tensor representing the bubble. Thus the integrals 
become more complicated. 

First we consider the 24 diagrams with a single internal gluon propagator, an example 
of which is given in figure 0(b). We write 



(C)na 2 = ~ / ix^P^ / dx x / '- Xl dx 2 / dLips[k - k u k 2 ]\M\ 2 C , (3.7) 

2! J k z i07T JO Jl-xi-e J 



1 rdk 2 C A C F 

y~ i ax i 

'l-xi- 

where the total amplitude contribution from these diagrams evaluated using the pinch 
technique is 

\M\ 2 



\C A C F N C ) eV ^ • (3.8) 

We again use the decomposition (|2.11 ) of C into primary, secondary and mixed terms, 
choosing to write 

<c«>»» = J c«(*7Q 2 ) ■ (3-9) 
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The details of the pinch technique require (see appendix [A]) that 
|M| 2 = {{h-hyg^ - 2k a g^ + 2k /3 g~< a } Tr 



1 1 UL 1 

Vila , : 7/3 7 , u TP27n 2 i // 7 ^ + ' 
Pi + f 1 Pi + f Pi + f 



(3.10) 

where the ellipsis represents 23 similar terms and the factor in curly braces arises from 
a modified triple-gluon vertex. (The contributions from the remainder of the triple-gluon 
vertex have been assigned in the pinch technique to other diagrams.) 

It is convenient to perform the integral over k\ and k 2 in the rest frame of k, in which 
we may write: 



k = Vfc 2 (l, 0,0,0) (3.11) 

ki = ^Vk?(l, sin cos 0, sin sin 0, cos 0) (3-12) 

k 2 = — sin 6* cos 0, — sin 6* sin 0, — cos^) (3.13) 

q = (g°,0,0,g 3 ) (3.14) 

Pi = (p°i,0,pl,pl) (3.15) 



and 



dLipsffc -> fci, k 2 ] = sin 6d6d<j)/(327r 2 ) . (3.16) 

It is in performing these integrals over 6 and <fi that the collinear divergence appears: in 
order to regulate it we substitute 



1 fc+fc h+h fa+h 



i>\ + h (pi + h) 2 (pi + ki) 2 - m 2 2pi ■ ki 



(3.17) 

v\=m? 



and similarly for the other terms. The integral then converges, with the divergences ap- 
pearing as terms in log A. 

In order to integrate over x\ and x 2 , it proves convenient to make use of the substitutions 

u = x\ — x 2 (3.18) 
e 

v H — = 2 — xi — x 2 = x 3 . (3.19) 

v 

After performing the integral over u we find we need to evaluate an integral over v, 

l dvf(v,e), (3.20) 

for some complicated function / involving logarithms and dilogarithms of rational functions 
of v and e. A direct expansion of the integrand in powers of e is not valid, but may be done 
by first dividing up the integration range 

f 1 dv f(v, e) = v r e f 1/K dyf^y, e) + [] dv f(v, e) , (3.21) 

where k is a small parameter, y/e < k C 1, on which the final result does not depend. (In 
practice we perform the calculation in the limit of small k.) 

The second term may now be evaluated by expanding in powers of e. Non-integer 
powers of e come only from the lower limit and always appear in the combination y/e/n. 
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In order to avoid having to keep track of terms e m {n/y/e) n , where m and n are large, we 
finally expand the result in k and neglect all terms that vanish as k — *■ 0. We obtain 



f 1 dv f(v, e) = a(log e) + ^6(log k) + (9(e) , (3.22) 



where a and b are polynomials. 

The first term of ( |3.21| ) may be evaluated by expanding the integrand directly in powers 
of e. This gives 



1/k 

dyf(Vey,e) 



6 (log k) 

■ + c 



K 



+ 0(e) , (3.23) 



where the constant c does not depend on e or k. Adding this to the result ( p.22|) gives the 
total contribution: 

/ dvf(v,e) =a(loge) + cV^ + 0(e) . (3.24) 

The individual components C^ 2 ( e ) are presented in the appendix; the total contribution 
from the 24 diagrams with a single internal gluon propagator was then found to be: 

C na2 (e) = -(4tt 2 - 33) log(eA) - — + ^C(3) - ( — log A + + 0(e) 

(3.25) 

where v ~ —89.556 is defined in equation (p.9| ). 
Note that: 

(i) The most divergent log 3 and log 2 terms found in the individual components have can- 
celled in the sum of contributions. There however remain log e and log A divergences, 
which will cancel if we sum all real and virtual diagrams. The leading divergence is 
again proportional to (47r 2 — 33), the leading-order value of (C). 

(ii) The coefficient of y/e in fl3.25p is simply the quantity c in equation ( |3.23j ). Thus an 
easy way to extract the coefficient of the 1/Q correction without performing the full 
calculation is to evaluate the integral 

= [ lK dy\imf(^fey,e) (3.26) 

Jl e^O 

1/k ( y 2 _ y dy r(y-l/y) r l f2n 



Ji 2y 2 J~(v~i/v) J-i Jo 27T e^o v 1 



704 Ji 2y 2 J-b-i/v) 

which is obtained by combining the definitions (|3.7|) and ( p.9|) with the changes of 
variable ( |3.16| ), ( |3.18[ ) and ( |3.19|) , and by making the substitutions v = -fky and 
u = \fez. The coefficient of the 1/Q correction is obtained by subtracting off the 
pole at k — 0, just as in ( |3.23|) . The integrand thus becomes identical with that 
obtained using soft matrix elements, as in [15|, [TJ|, but the phase-space boundary is 
more carefully defined in terms of the small arbitrary quantity k. 

The final diagrams are those such as figure 0(c), diagrams which have no internal gluon 
propagator. There are 36 such diagrams, 18 of which have colour factor C\ and 18 of which 
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have Cp — ^CfCa- Let us choose to write the contributions to the mean value of the C 
parameter from these diagrams respectively as 

c F rdk 2 nc F 



(r (i)\ C F \ dk LLC F 2 r (i) ( h 2 jn 1 



= %\ ^ 11(2 ^" CA) a^(fc a /Q a ) • (3-27) 

The integration is performed in the same manner as above. The individual components 
are listed in appendix |TJ, and we find evaluate the total characteristic functions to be: 

CUe) = -Vi(^logA + ^ + ^+0(c) 
CUe) = -^(4vr 2 -33)logA + ^-^C(3) 

+Vi(^logA+^ + ^+0(e) (3.28) 

Again the leading divergence is as expected proportional to (Air 2 — 33). 

Now that all the diagrams have been calculated, we are in a position to write down the 
total real contribution to the 1/Q correction to (C). 



3.3 Total real contribution 

The total real contribution to the mean value of the C-parameter, expanded in powers 
of e, is, by simple compilation of the results above, 

w In J k 2 I An n \V J p \ 704 3 128 3 R f J J 

HCU 2 /l2 2 7 2 8 64 1 27^-3^1 



" F { (An 2 - 33) log A - 54 + 48((3) } + 0(e) 



TT 



fc7Q 2 , (3.29) 



where the first line comes from the insertions into the gluon propagator, (|2.16|) and (|3.5|) , 
and the second and third lines from the remaining diagrams ( |3.25| ) and (|3.28| ). 

The factor (— /3 Q / An)\a s (— k 2 )\ 2 is naturally interpreted as the spectral function p s (k 2 ), 
generated by making renormalon insertions into the gluon propagators. It is not clear that 
the same assignment can be made with the remaining contributions, since there is not a 
simple diagrammatic way of generating the running in diagrams without propagators in 
which to insert renormalons. If naive non-abelianisation gives the terms in the first line of 
( |3.29| ), the terms on the second and third lines are extra. Nevertheless, the Milan factor is 
calculated on the assumption that this assignment is valid, and thus we obtain for the real 
part of the Milan factor 



4 15n 2 ' 

Meal = i-^CU r"""" +^-1 -^T R n f - 




704 J 3 J 128 
l0n 2 T R n f 



64 \\C A -AT R n s 



(3.30) 
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4 Discussion 



In the preceding sections we have calculated the mean value of the C-parameter in e + e~ 
annihilation from 2-loop real emission diagrams as a function of the invariant mass k 2 = eQ 2 
of the two secondary partons. The coefficient of the y/e term in the series expansion is then 
assumed to determine the magnitude of the 1/Q power correction. 

There are a number of points to note here. The first is that virtual contributions have 



not been calculated exactly. We include them here in the same manner as in []To], |I6j| , using 
the soft approximation: 

M ,„ = _ W-W)c* 



Details of how this number is extracted from ]15| are found in appendix 0. Thus we predict 



the total 1/Q correction to the C-parameter to be 



5(C) = 12M^ , (4.2) 



where the 2-loop enhancement factor is: 



3 (128vr + 1287rlog2 - 35n 2 )C A - lOir 2 T R n F 
~ 64 11C A - 4T R n F ^ ' ' 



Putting rif = 3 yields M. = |vr(l + log 2) — -^tt 2 m 1.49, in agreement with the existing 



numerical calculations [IB], M. The use of the soft approximation in calculating such 
enhancement factors is therefore shown to be valid. 

It may not at first sight be totally clear exactly how the virtual diagrams, which do not 
have any phase-space restriction, could contribute to the 1/Q correction in the full two- 
loop approach adopted here. However a term proportional to \fjj? /Q, where fi 2 is a dummy 
variable in the loop integral, might well appear even without a phase-space restriction if 
the integrand happens to be of the correct form. 

It is also not clear exactly how the spectral function p s {k 2 ) is generated diagrammatically 
- in the large nf limit this is achieved by quark-loop renormalon insertions into the gluon 
propagator, which can be extended to gluon-loop insertions using the pinch technique. 
This reconstructs the correct value for ft , and is a diagrammatic expression of naive non- 
abelianisation. But there are left-over contributions to which this does not apply. Whether 
or not these can truly be said to contribute on the same basis as the renormalon-type 
contributions is an open question. 
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A Appendix — diagram decomposition 



We use the pinch technique [[?[]- [|I0| to reallocate certain gauge-dependent terms between 
diagrams. We first write the diagrams in the simplest gauge for the process, a naive 
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Feynman gauge with external ghost fields to cancel the unphysical degrees of freedom, 
and then apply pinch technique manipulations and re-assignments. Of course by gauge 
independence the same results may be obtained with any gauge choice (for a detailed 



analysis, see 11 



The diagrams in figure ||(c) are (except for the colour factor) the same as those that 
appear in QED. With this gauge choice these diagrams become simply: 



(C F N C ) eV Tr 



f>H c 



-r 



+ ■ 



(a.i; 



where the ellipsis represents the remaining 35 diagrams. We calculate these diagrams just 
as they are, without further modification. 

The diagrams containing a single triple gluon vertex and gluon propagator, such as in 
figure 0(b), are given in this gauge by: 



(\C A C F N C ) eVp {(h-hY'g^ - (k + k 2 ) a g^ + (Ar+fcO V°} 



Tr 



k 2 



X 



"7/3 " 



;7 7 + ■ ■ 



% + h Ip h + r h + 

Here we split the triple gluon vertex into two pieces: 

{(h-h^g" 13 - {k + k 2 ) a g^ + (fc+fci)V a ) = {(h-h) 1 'g aP - 2k a g^ + 2fcV a } 



(A.2) 



{Kg 



Pi 



}■ 



(A.3) 



The contribution from the first term on the right hand side is that used in the calculation of 
the diagrams, equation ( |3.10| ). The contributions from the second term (such contributions 
are called the 'pinched parts' of the diagrams) are manipulated using identities such as: 



A o 



1 



-7/3 ; 



1 



to bring it finally into the form 

4 (C A C F N C ) eV (k 2 g al3 - k a kp) Tr 



i>i + 



77^ + 



(A.4) 



(A.5) 



The Dirac trace appearing here is identical to that found in the calculation of diagrams 
H|(a): it contains 4 terms in total. This contribution will be added in due course to the 
conventional loop insertion diagram. 

These loop insertion diagrams, figure @(a), are given by 



(C A C F N c )eY^ {(k 2 -hyg^ - (fc+A; 2 )V p + (k+hfg^} 



x 



(A.6) 



{(k 2 -k 1 ) a g a/3 - {k+k 2 ) a gp„ + (k + k^pg^ Tr 
which reduces to 



1>i + $ 



tl+k 1 



7 CT + 



{C A C F N C ) eVrj {A(k 2 g aj3 - k a kp) + 2(k la - k 2a )(k w - k w ) - (k la k 2f3 + k 2a k lf3 )} x 



Tr 



(A.7) 



12 



The only remaining contribution is from the cut ghost loop, which may be written 
1 



(C A C F N C ) e g — {k la k 2/3 + k 2a hp} Tr 



and thus the total contribution to the cut loop insertion diagram, using this prescription 
for rearrangement of gauge-dependent terms, is just the sum of ( A.5 ), ( A.7 ) and (|A.8 ): 



(C A C F N C ) {8(k 2 g a p - k a k p ) + 2{k la - k 2a ){k w - k w )) x 



Tr 



Pl + P Pl + f 



(A.9) 



It is this combination of terms that is represented by the diagram 0(a) and the calculation 
corresponding to it. 



B Appendix — individual contributions to (C) 

Here are presented the individual contributions to (C), denoted C^(e) in the text: 
(i) In the massive gluon scheme, we find: 



Cmg( e ) 



31og 2 e + 91oge-vr 2 + 15 + 0(e) 

223 



3 log 2 e - 171oge + vr 2 - ==p + 0(e) 
65 



" " ' " 6 

„ 65 

.loge + 47r 
= 0(e), 

which give the total 

C mg (e) = 4tt 2 - 33 - 12^ + 0(e) . 
(ii) For gluon splitting into a secondary qq pair, we obtain: 



C. 



c. 



ab ^ 

(m) . 
ab I 



3 log 2 e + 9 loge - Tr 2 + 15 + 0(e) 

223 

-3 log 2 e - 17 log e + Tr 2 + 0(e 

6 

, „ 65 45vr 3 n m , , 
8 log e + 4vr 2 - — - — - ^ + 0(e) 
6 32 

0(e), 



with total 



Cab(e) 



4^-33-^^+0^). 



T3.1) 



(B.2) 



(B.3) 
(B.4) 
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(iii) From the cut loop diagrams such as in figure |](a), we get: 

Cffi(e) = 31og 2 e + 91oge-vr 2 + 15 + 0(e) 
C&(e) = -31og 2 e-171og e + 7r 2 -^P + 0(e) 

&{e) = 81og6 + 47r 2 -|-^^+0(6) 

C { :ke) = 0(e), (B.5) 



with total 



261 7T 3 

C na i(6) = 4tt 2 - 33 - -=rV~e + 0(e) . (B.6) 

176 



(iv) The diagrams such as in figure |2|(b) yield 

( e ) = yt log3 e + IT logA log2 e + IT log2 Aloge + YT log2 e 

108 , A , 36 , 2 a 30tt 2 - 261 , 12tt 2 -171 . 
+Yj-logAloge+ — log A — loge — logA 

39tt 2 477 270 A/ x 12tt 2 , 24, , 144 T ... ^. . 

+ C 3 + log 2 bg 3 2 + L 3 i + 0(e) 

11 11 ll sv; 11 s 11 6 11 SK2J yj 

~(n\ / s 12 , q 36 , . , o 18 , n . , 102 , r, 

Ci (e) = - log 3 e - — log A log 2 e - — log 2 A log e - — log 2 e 

204, 54, 2 A 30tt 2 - 507, 12tt 2 - 409 , . 

— rr g ge ~TT g — n — lo s e + — ^ — lo § A 

9 1042 270,. « 12tt 2 , 24, , 144 ... 

+7tt 2 - — + — C(3) - — log 2 + - log 3 2 - —Lis® + 0(e) 

j ml/ , 48, 2 96, . 18 2 a 48tt 2 -150, 
ffi(e) = J! log e+ IT logAloge + — log 2 A + loge 

48tt 2 - 158, . 38tt 2 785 1440 >/oN 

+ log A + C(3) 

11 6 11 11 11 sv ; 

^/367r 2 , . 144tt 2 \ _ . 

~^[ir g + ^t~ + z/ ) + (e) 

C&(e) = 0(e), (B.7) 
which have total 

„ / x 12 / 2 m / ax 1350 1440 A/ , 
C na2 (e) = -(47r 2 -33)log(eA)- Tr + ^ r C(3) 

n /367T 2 , . 1447T 2 \ , nl , .„ , 

-^f— logA + ^^ + z/j +0(e) . (B.8) 

The number z/ « —89.556 is given by the expression 

288tt r&r/e . 4327T . 967T 2 lQg2 
f = — — / log sin -dd — / log sin . (B.9) 

11 A/6 2 11 A/3 2 11 
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(v) The diagrams such as that of figure |2|(c) which have colour factor C\ give the contri- 
butions: 



Mm). 
u na3 \ 



e. 



na3 ' 



18 , 2 a i 18 1 2 a 9 , 2471 " 2 - 9 

— log A log e log AH log e H 

n » s n 6 T 22 s 22 

^ 108 
11 22 11 sv ; w 



log A 



18 1 2 A , 48 1 2 A 

— log 2 A log e H log 2 A 



73 , 24vr 2 - 73 , . 
— log e log A 

22 6 22 & 



7 



30, n 32 , 32 , . 10vr 2 

— log AH log e log A H 1 

n s t 1X s 11 s 11 11 

^/36tt 2 . 36tt 3 \ , 
--v/el — \og^ + - rr + v) + 0(e) 



11 



11 



0(e), 



(B.10) 



which have total 



C n a3( e ) 



'36vr 2 , . 36tt 3 \ _ . 
- r I —r^ log A + + v]+ 0(e) . 



11 



11 



(B.11) 



(vi) And finally those diagrams with colour factor C\ — \CaCf provide: 



^na4( e ) 



^nai ( e ) 



^na4( e ) 



36 2 I 8 i 2 a i 

-yj log A log e-— log A loge 



108 , A 1 36 1 2 A 

— log A log e — — log A 



36tt' 



11 



135 , 12vr 2 
log e + - 



171 , . 48tt 2 711 252 A/ , , nl . 

log A + + C 3 + 0(e 

11 6 11 22 ll su u 



36 , a , 2 18. 2 a i 204 i A i 54 2 a 

— logAlog e + — log Aloge+ — logAloge + — log A 
36vr 2 -219, 12vr 2 - 409 , . 86vr 2 1487 252 



11 



loge 



11 



log A 



96 , a i 18 i 2 a 

■— logAlog e - — log A 

38vr 2 260 576 „ x 

— + TT-TT« 3 » + - 



11 

84 , 48vr 2 
loge - - 



22 
158 



11 



C(3) + 0{e) 



0(e) 



11 

^1 — logA 



11 

36tt 3 

11 



log A 



v 



0(e) 



(B.12) 



with total 



„ . . 396 - 48vr 2 , . 648 - 576C(3) 
C na4 (e) = - logA + - 



11 

r(^\ a 



11 

36tt 3 
11 



0(e) 



(B.13) 



C Appendix — extraction of virtual contribution 



In order to compare our real-emission result ( |3.30| ) with the Milan factor we have to 



note that the latter as derived in [[Lq, 16] includes a contribution from virtual diagrams as 
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well as real diagrams. This is a part of the inclusive piece r inc i = c^Ca/Po, where from 



equation (3.8) of |TJ| we have 
2x 



,(0 



dx 



— — - (y/l + x 2 - x) log 

1 + X 1 v ' 



x 2 (l + x 2 ) 



8-41og2-27r(l -log 2) . (C.l) 



Here the contribution involving yl + x 2 is the contribution from real emission diagrams, 
while the virtual contribution is represented by the — x term. The combined integral con- 
verges at large x. 

In the full calculation (i.e. without using the soft approximation) the real and virtual 
contributions represented here arise from integrals of the form: 



A^re/vi = / d£ /i rc /vi(£,e) 

J 



(C.2) 



where the kinematical integration variable is £ — k±/Q. Using ( |3.1| ) and the argument of 
| HI we have in the soft (small £) region that 



M£,e) - 

M£> e ) 



2t 



£ 2 + e log 



Olog 2 



e 2 + e 



(C.3) 



It is therefore clearly not permissible to evaluate the integrals /C re / V i by a naive expansion 
of the integrand in powers of e. Instead we divide up the integration region according to 



rl/K rl 

£re/vi = \/e/ dxh ve/vi (y/ex } e) + d£h ie / vi (£,e 

JO Jy/e/K 



(C.4) 



where the total cannot depend on k. 

Neglecting terms that vanish as k —>■ 0, the first term on the right-hand side gives 



V« , 2x 



dx vl + x 2 log 

1 + x 2 

i/k 2x 2 

dx— — g log 
1 + x z 



x 2 (l + x 2 ) 



x 2 (l + x 2 ) 



(l-log«) + 8-41og2 

K 



K 



:i-log«)-27r(l-log2) , (C.5) 



while in the second the \fi terms can come only from the lower limit, and thus appear 
as \/e/ k. The 1/k terms in ( P.5|) are therefore exactly cancelled, leaving the total y/e 
contributions from real and virtual parts respectively as 



-Teal 



8-41og2 



-2tt(1 - log 2) . 



(C.6) 
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